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Abstract 

We study the quantum mechanics of a generahzed version of the baker's 
map. We show that the Ruelle resonances (which govern the approach 
to ergodicity of classical distributions on phase space) also appear in 
the quantum correlation functions of observables at different times, and 
hence control the statistical variance of matrix elements of observables 
(in the basis of eigenstates of the quantum time evolution operator). We 
illustrate this with numerical results. 
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I. INTRODUCTION 



The baker's map and its generalizations [|T],0 provide prototypical models of chaotic sys- 
tems, and thus their properties, both classical and quantum, have been extensively studied. 
The issue we examine here is the structure of the Ruelle resonances (which control the 
approach to ergodicity of smooth distributions on the classical phase space), and the role 
these resonances play in the properties of the corresponding quantum map. 

This question has been partly addressed in earlier work. In ref. [Q, the Ruelle resonances 
of the original baker's map were computed, and the classical dynamics of phase-space distri- 
butions were thoroughly investigated. In ref. |^ , the statistical distribution of the eigenvalues 
of the quantum time-evolution operator U for a generic chaotic system was studied, and it 
was found that the Ruelle resonances control the deviations of this distribution from the 
Wigner-Dyson form that is predicted by the random- matrix analogy 0. In ref. [0, certain 
matrix elements of the resolvent of U were studied (for the original baker's map), and used 
to identify particular Ruelle resonances. Our work here is closely related, though we require 
less numerical effort to see the resonances in the quantum time-correlation functions that 
we examine. We also uncover the role of the resonances in the statistical variance of matrix 
elements of observables (in the eigenbasis of U) . 

We review the needed information about the classical map in Sec. y, and about the 
quantum map in Sec. |IT|. In Sec. |^we introduce the quantum correlation functions that we 
study numerically in Sec. |V|. Sec. [V^ contains our conclusions. In an Appendix, we calculate 
the values of the Ruelle resonances by computing the associated spectral determinant or 
dynamical zeta function. 



II. PROPERTIES OF THE CLASSICAL MAP 

Classically, the generalized baker's map Ai: x ^ x' = Ai{x) is defined on a unit 
square (interpreted physically as a two-dimensional phase space) x = {q,p), < q < 1, 
< p < 1. To implement the map, this square is first cut into s vertical strips with widths 
Wi, . . . , Wg] J2l=i Wa = 1- Let Ca = J2b=i '^b bc the left edge of the a**^ strip (with ci = 
and Cs+i = 1). Each strip is then stretched horizontally and compressed vertically to unit 
width and height Wa, and the strips are stacked vertically to assemble a new unit square. 
We therefore have 

= 2^\ , ^(ea <q< ea+i), (1) 

a=l \ WaP +ea J 

where we have introduced a generalized step function: 

Note that if q is in the range of the a^^ strip (ca < g < Ca+i), then so is p' (cq <p'< Ca+i). 
This point is important in the construction of the quantum map. 

It is convenient to introduce a Dirac notation: let |a;) be a phase-space eigenstate with 
normalization {x\x') = 6{x — x') = S{q — q')6{p — p'). (To avoid confusion we use parentheses 
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for classical states and angle brackets for quantum states.) We define the Perron- Frobenius 
time evolution operator U via 

{x\U\x') = 5{x-M{x)). (3) 

If we consider the domain of U to be the space of Lebesgue square-integrable functions on 
the unit square, then U is unitary for any area-preserving map M.^ and so its spectrum lies 
on the unit circle. However, if the domain of U is restricted to a suitable class of smooth 
(infinitely differentiable) functions, then U is effectively truncated, and it is no longer unitary 
in the reduced space. If this truncated lA is put into Jordan (upper-triangular) form, the 
diagonal entries are inside the unit circle; these diagonal entries are the Ruelle resonances 
of the map. Their values control the decay of initially smooth distributions on phase space 
to the ergodic (uniform) distribution. 

Hasegawa and Saphir have shown that the Ruelle resonances of the original baker's 
map (s = 2, wi = W2 = \) are given by ui = with degeneracy di = £ -|- 1 for 
i = 0,1,2,.... Below, we show that for the generalized baker's map, ui = J2l=iwl'^^, 
with the same degeneracy. In general, the values of the Ruelle resonances have no direct 
connection with the Lyapunov exponents that govern the short-time instabilities of individ- 
ual trajectories. For example, the positive Lyapunov exponent for a generic (nonperiodic) 
trajectory is A = — J2l=i Wa Inwa- 

Consider the time evolution of an initial smooth distribution on phase space, A{x) = 
{x\A). We have 

{x\U'\A) = J dx' {x\U'\x'){x'\A) 

= J dx' 6{x - M'{x'))A{x') 

= A{M-\x)). (4) 

Here denotes t iterations of A4, and A^~* denotes t iterations of the inverse map A4~^] 
the last equality in eq. (§) follows from the change of variable x' = (?/'), which has unit 
jacobian if Ai is area preserving. 

Let us now consider time correlation functions of the form {B\U'^\A), where A{x) = {x\A) 
and B{x) = {x\B) are both smooth functions, and t is an integer (positive or negative). We 
have 

iB\U'\A) = Jdx {B\x){x\U'\A) 
= J dxB{x)A{M-\x)) 

= J dxB{M'{x))A{x). (5) 
For the generalized baker's maps, these correlation functions can be written, for t > 0, as 

oo 

{B\U'\A) = J2ce{t)ui (6) 

where Ci{t) is a polynomial in t of maximal order di — 1 = i. The coefficient Ci{t) depends 
on the choice of A and B, but the Ruelle resonance ue does not. We always have uq = 1 
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and Im^I < 1 for i > 1, and so initial correlations decay as time evolves. If we consider the 
case A = B, then eq. implies that we have time symmetry, and so 

oo 

{A\U'\A) = T.cmnf, (7) 

e=o 

where C£(|t|) is again a polynomial of maximal order i in \t\. 

It is useful to introduce a set of basis polynomials (actually shifted and rescaled Legendre 
polynomials) M 



that are orthonormal on the unit interval, Jq dq Pe'{q)Pe{q) = See. We define a set of classical 
states \im) on phase space via 

{x\im) = Pe{q)Pm{p). (9) 
According to eq. (HI), the action of the Perron- Frobenius operator on these states is given by 



[X 



U\em) = J2 Pi{Waq + ea)Pm{iP - ea)/Wa) e{ea <P< Ca+l). (lO) 



a=l 



We wish to find the matrix elements of U in the \im) basis, 

{e'm'\U\im) = J dxPe'{q)Pm'{p){x\U\em). (11) 

In order to do so, we must evaluate integrals of the form Ip£ = Jq dq P£'{q)Pe{wq + e), where 
< w < 1 and < e < 1 — w are constants. (The p integrals can be put into this form by a 
change of variable that results in an extra factor of w.) Using eq. (|) for P^'(g) and repeatedly 
integrating by parts, we see that If/e = if £' > £, and it is similarly straightforward to show 
that III = w^. Thus the matrix {i'm'\ U\im) is in Jordan form (upper-triangular with respect 
to the i indices, and lower-triangular with respect to the m indices). The diagonal elements 
are the Ruelle resonances, given by 

{im\ U\im) = ^ = ue+^, (12) 

a=l 



in agreement with eq. (|A9|) . We also see that ue has degeneracy d^ = i + 1, since there are 



+ 1 different \im) states that result in the same value of {im\U\im). 
Eq. (12) also holds if we replace U with U^^ on the left-hand side. Furthermore, because 



of the Jordan form of U in this basis, we have 

{im\U'\im) = (13) 

for any integer t (positive or negative). Thus the time-evolution of one of the \im) states 
picks out a particular Ruelle resonance. 
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III. PROPERTIES OF THE QUANTUM MAP 



To quantize the generalized baker's map, we must discretize q and p let qj = {j—vi)/N 
and pj = {j — 1/2) /N, where j = 1,. . . ,N] the integer N plays the role of the inverse of 
Planck's constant h, while < < 1 are parameters of the discretization that should be 
irrelevant in the classical limit. Also, we require that the strip widths Wa be integers divided 
by A^. We then define corresponding quantum states \qj) and \pj) with the properties 

{qM = {Pj\pk) = 5jk, (14) 



{Pj\qk) = N-^/^exp{-2mpjqk/N) = {F, 



N)jk- 



(15) 



We now wish to specify the unitary time evolution operator U for a single iteration of the 
map. From eq. (|l|), we expect that, heuristically, U\qk) ~ \{qk — ea)/wa) when is in the 
a*^ strip. This must be reconciled, however, with the evolution of the momentum p. Recall 
that if the initial coordinate q is in the a^^ strip, then so is the final momentum p'. This 
motivates a specification of in a mixed p-q basis. 



{Pj\U\qk) = J2^a ^^"^e'^'iPjUQk - ea)/wa) 9{ea < qk < Ca+l) 0{ea < Pj < Ca 



(16) 



a=l 



The prefactors of w~^^'^ are needed for unitarity, while the phase angles Xa can be arbitrary 
linear functions of qj and pk', these phases should be irrelevant in the classical limit, since 
they can be absorbed by shifts in the origins of q and p in the definition of the inner product 
(p!5|). It is then convenient (and by now traditional) to choose these phases so that U is 
given in the coordinate basis by 



{Qj\U\qk) 











Nw2 







where the 0a's are arbitrary constant angles. 
Since U is unitary, its eigenvalues are phases, 

U\a) = e-'^-\a), 








ik 



(17) 



(18) 



and the corresponding eigenstates are orthonormal, {a\ 



IV. QUANTUM TIME CORRELATION FUNCTIONS 

Consider a quantum operator 0{q,p) that is a smooth function of the position and 
momentum operators q and p, and that does not depend explicitly on h = 1/2'kN. Then 
Shnirelman's theorem states that the diagonal matrix element where \a) is an 

eigenstate of U, tends to the phase-space average / dxO{x) as h 0. (For a proof in 



5 



the case of the baker's map, and a complete set of references, see 0.) For hamiltonian 
systems, this result was generalized by Feingold and Peres JlOl and by Wilkinson (for 
rigorous results see |]12[) to the case O = BU^AU~^, where A{q,p) and B{q,p) are smooth, 
A^-independent functions. In the Heisenberg picture, U^A{x)U~^ = A{xt), where Xt denotes 
the time-evolved position and momentum operators. If, however, we use classical evolution 
instead, then according to eq. @ Xt A^^*(s). This substitution should be valid in the 
N —>■ oo limit, and so 

{a\BU'AU-'\a) = J dx B{x)A{M-'ix)) + 0{N-^/^), (19) 



where the 0{N ^/^) corrections can be expressed as a sum over periodic orbits ||TT] , p!2 



If we insert a complete set of U eigenstates on the left-hand side of eq. ([I9|) , we have 
{a\BU'AU-'\a) = J2{a\BU'\(3){(3\AU-'\a) 

= Y.B^pA^^e^^'--'^^\ (20) 

where B^p = {a\B\(3), etc. Thus we get the Feingold- Peres- Wilkinson formula (applied to a 
quantum map rather than a hamiltonian system), 

J2B^pAf^^e'^'--'^^' = I dxB{x)A{M-\x)) + 0{N-^'^). (21) 

Now, according to eq. (^, the leading term on the right-hand side of eq. (|2l|) is given by 
the right-hand side of eq. (^. Thus we have, for t > 0, 

oo 

^ B^pAp^ e'^'--'^^' = ^ Q(t)e-^^* + OiN^'/') (22) 

/3 £=0 

where we have defined e"''''^ = U£. This is our main result, which we will investigate numeri- 
cally in the next section. 

Also, since the leading term on the right-hand side of eq. (|2^ ) is independent of \a), we 
can improve the accuracy (by a factor of N^^^"^) by averaging over a, 

-1 oo 

- Y: BapApa e'^'--'^^' = Y: Q(t)e-^^* + 0{N-'). (23) 

Q/3 e=o 

Note also that the left-hand side of eq. (|23|) can be expressed as the quantum trace 
{l/N)TrBU^AU-\ 

Let us consider the special case A = B; we then have 

-1 oo 

^$:|A„^pe^(^'^-^^)* = 5:Q(|t|)e-^^l*l+0(Ar-i). (24) 

We see that the statistical variance of the matrix elements A^js is controlled by the Ruelle 
resonances. For example, let us choose A{x) = Pe{q)] then a single Ruelle resonance appears 
on the right-hand side of eq. (plD, and furthermore Q(|t|) = 1 in this case. Thus we have 
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1 Y: mq)ap\' e^^'--'^^' = e-^^l*l + 0{N-'). (25) 

We next multiply both sides by e"*'^*e"^l*l, where < e <^ min(l,7£) provides a cutoff, and 
sum over t; we get 

N 5 sin2((^. -Op- co)/2) + (5/2)2 = ,1^2(^/2) + srf(7./2) + ^^^''^^ ^^^^ 

If we now take the formal limit of e — > on the left-hand side, the e-dependent factor becomes 
27T6{da—6i3—uj). This shows us that the statistical distribution of the matrix elements Pi{q)ai3 
with fixed u = 6a — df3 is a Lorentzian in sin(ti;/2) whose width is 2 sinh(7^/2). However, this 
formal limit must be interpreted with care, and so in our numerical work we concentrate on 
eq. (HD. 



V. NUMERICAL RESULTS 

We consider the example of a two-strip map with wi close to 1/3 and W2 close to 2/3. 
We then numerically evaluate 

CeAt) = {a\P,{q)U'P,{q)U-'\a) = ^ mqU\' e'^'"-'^^' (27) 

and its average 

Ceit) ^j-TTPeiq)U'Peiq)U-' = 1 ^ |P,(g),^|2 e^C^-^/')* (28) 

a/3 

for different values of £ and A^. We note that the numerical problem of computing Ce{t) 
is considerably simpler than that of computing Ce^t); to compute the latter, we must 
diagonalize U, whereas the former can be done by taking a trace of a simple matrix product 
in the position representation. We take the irrelevant parameters to be z/i = 1/2 = | and 
01 = 02 = 0. We compare these results with the classical expectation e""^*^'*' in fig. (1) for 
= 331, the highest value we used; good agreement can be seen. 

Because of the 0{N~^) corrections in eq. (^4]), agreement between Ci{t) and e~^'^^^ breaks 
down at later times. For a given choice of i, this time is estimated as 

log N 

tiog = • (29) 

Physically, for t > tiog, classical evolution produces phase-space structures with areas less 
than h = 1/N, and so quantum evolution is necessarily different. We can define a breakdown 
time in practice as the first value of t for which C^(t) > Cg{t — 1). The results are compared 
with eq. ( pOf ) in fig. (2) with good agreement. 

It is also interesting to look at off-diagonal correlation functions. Let us consider the 
case A{x) = (x|10) = Pi{q) = ^3(1 - 2q) and B{x) = {x\01) = Pi{p) = ^3(1 - 2p). 
Classically, we can straightforwardly evaluate the single-step correlation (01| W|10); the result 
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is 1 — J2l=iwl = 1 — U2- Also, (01|10) = by orthogonality. Since only the resonance ui 
can be involved in the time correlation of these distributions, and since the coefficient Ci{t) 
can be at most linear in t (because Ui is two-fold degenerate), we must have 



(01|W*|10) = {1-U2)tu\ 

for t > 0. (For any t < 0, this correlation vanishes.) 

Quantum mechanically, it is also easy to evaluate the t = 1 correlation function, 



(30) 



1 

N 



TrPi(p)f/Pi(g)[/-i = -Y^p,{p^){p^\U\q,)P,{q,){q,\U-'\p,) 



N 



jk 



jk 



1 Wn. Jen. Je„, 



a=l 
s 

E 



1 



[V3(ea+i 



a=l 

i-E:=x+o(iv-i^ 



0{N 



-1^ 



(31) 



In the third line, we used eq. (|T6[) and took the N ^ oo limit. To get the last line, we 
repeatedly used Yfa=i{.^^a+i ~ ^a) — ^ exponent z/ > 0. Eq. (^) gives a simple 

example of how quantum expressions go over to their classical limits. 

In fig. (3), we show the quantum and classical expressions for this correlation function, for 
our two-strip example. To (slightly) improve the agreement, we also performed an ensemble 
average over ten random choices of the four irrelevant parameters. (This can also be done 
for the diagonal correlation functions, with a similar slight improvement.) 



VI. CONCLUSIONS 

We have examined, both analytically and numerically, quantum time correlation func- 
tions of the form Tt BU^AU~^ in a generalized baker's map, where U is the quantum time 
evolution operator, and A and B are smooth, h independent, functions of the position and 
momentum operators. In the ^ — >■ limit, these quantum correlation functions are governed 
by the Ruelle resonances that govern the approach to ergodicity of the classical distributions 
A and B on phase space. We see this very clearly in our numerical results as long as the 
correlation function itself remains larger than the predicted 0{h) corrections. These results 
also imply that the statistical variance of the matrix elements of a smooth observable (in 
the eigenbasis of U) is controlled by the Ruelle resonances. 



APPENDIX A: RUELLE RESONANCES FOR THE GENERALIZED BAKER'S 

MAP 

The Ruelle resonances ue = e~^^ of the generalized baker's map are given by the zeros of 
the spectral determinant 
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d{z) = det{l- z-^U). 



(Al) 



We compute d{z) following the analysis of Hasegawa and Saphir for the original baker's 
map. We write 



d{z) =exp[Trln(l-2-iW)] = exp( - ^ — Tr^^' 



oo ^-t 



t=l 



(A2) 



which is valid for \z\ > 1. The trace we need is given by 



TrW* 



dx S{x - M\x)) = dqo dpo S{qt - go) ^(Pt - Po) 



(A3) 



where qt and pt are the values of q and p after t iterations of the map. To evaluate this trace, 
we must find the periodic orbits of period t. Each of these orbits is uniquely labelled by the 
set of strips in which we find qo,qi, . . . , qt-i] the total number of these orbits is s* (where s 
is the number of strips). The trace can then be evaluated as a sum over the periodic orbits, 



TrW* = ^ 



p.o. 



det 



' dxt 



dxo 

E|(A-^-l)(A-l) 



p.o. 



= EE(^ + l)A^^^ 

p-o. e=o 

where we have defined 

A = dpt/dpo = {dqt/dqo)~^ = w{qo) . . . w(gt_i). 



(A4) 



(A5) 



and w{qi) is the width of the strip in which we find g^; A is related to the positive Lyapunov 
exponent A for the orbit by A = e"'*'*. 

Let rii be the number of times q is in strip i during a particular periodic orbit; J2i=i = t- 
Then, for that orbit, A = w"^ . . . and there are t\/{ni\ . . . n^!) orbits with this value of 
A. Thus we have 



E---E 



^ ^ 'ni+...+ns,t 

ni=0 ns=0 '^l- • • • e=0 



(A6) 



We exchange the order of the summations to get 

oo i t j.| 

e=0 ni=0 ns=0 'n- ■ ■ ■ It's- 

oo 

= E(^ + i)«' + --- + ^f')*- 



(A7) 



Inserting eq. (|A7|) into eq. (|A2|) we get 
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/ oo -t oo ^ 

ci(z)=exp -E — , (A8) 

^ t=l 1=0 ' 

where we have defined 

= wf+i + ... + M;f+i. (A9) 
Again exchanging the order of the summations we find 

d{z) = exp ( f;(£ + 1) ln(l - z-^uS] , 

oo 

= \\{\-z-^u,)'^\ (AlO) 

If we now analytically continue to \z\ < 1, we find zeroes of order i + 1 at z = Ui for 
£ = 0,1,2, ... . This demonstrates that the Ruelle resonances are given by eq. (|A9|) , and 
that ue has degeneracy de = Note that uq = 1; the corresponding eigenfunction is the 

ergodic distribution, which is uniform over the unit square. 
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FIGURES 




FIG. 1. Points show the correlation function Ci{t) = (l / N) Ti Pi{q)U*Piiq)U'-* for 
£ = 1,2,3, with N = 331 and wi = 110/331; error bars show the root-mean-square range of 
Q,a(i) = {(^\Pe{Q)U^ Pe{Q)U~^\a) . SoHd hnes show the corresponding classical correlations u\. 
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FIG. 2. Points show the break time tiog, defined as the first value of t for which the correlation 
function Ce{t) increases, as a function of I = 1, 2, 3, for N = 21, 72, and 331. Solid lines show the 
the prediction of eq. ([29|) . 
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FIG. 3. Points show the ofF-diagonal correlation function {l/N)Tr P£{p)U^Pi{q)U~* for 
N = 331 and wi = 110/331, ensemble-averaged over ten random choices of the classically ir- 
relevant quantization parameters. The solid line is the classical expectation (1 — U2)tu'^^ ■ 
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